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SECTION 1 - INTRODUCTION ‘ ‘ ' i 

A dependable environmental data base is required by the agriculturalist in 
his crop management practices and problems. Much of this is provided directly. 
An important indirect channel is computerized crop management models such as 
those developed by Virginia Polytechnic Institute and State University (VPI5SU), 
aimed at optimizing irrigation procedures and controlling plant diseases and 
pests. These models require reliable estimates of precipitation on a day-to-day 
basis in near real time. Daily observations of precipitation reported on a 
network of rain gages provides accurate data for the immediate vicinity of each 
gage. If the distribution of gages were uniform and sufficiently dense (say 
several hundred gages over the state of Virginia) and if all measurements were 
made at the same time of day, and made immediately available to the computer, 
then interpolation between gages should provide a reliable estimate field for 
the entire state. In practice, meeting these criteria is both difficult and 
expensive. Furthermore, convective rain, the predominant type during the 
growing season, is usually sharply discontinuous, frequently invalidating 
straightforward interpolative methods. 

A number of investigators have developed and tested methods to estimate pre- 
cipitation using meteorological satellite data. See Barrett^ \ Chan^ , 
Follansbee and Oliver Follansbee^^\ Griffith et al.^^\ Martin and 
Scherer Scofield and Oliver and Woodley and Sax^^^. These techniques 
show varying skill in delineating precipitation patterns in the data-sparse 
areas between gages. Geostationary satellites, which provided complete cloud 
photograph coverage of the United States every half hour during daylight hours 
in the visible spectrum and around-the-clock in the IR, are particularly useful. 

An excellent technique based on geostationary satellite data developed by 

f81 

Scofield and Oliver , has been used successfully by the Synoptic Analysis 
Branch of the National Environmental Satellite Service (NESS) in support of 
the Quantitative Precipitation Branch (QPB) and the Weather Service Forecast 
Office (WSFO) of the National Weather Service (NWS) in estimating heavy 
precipitation and flood situations. Another such technique, developed jointly 
by the University of Wisconsin and the National Hurricane and Experimental 
Meteorology Laboratory, Miami, Florida, has been used in support of the GATE 
project of the Global Atmospheric Research Program. See Martin et al.*’^^^ 
and Griffith et al.'" . Both methods depend in large part on the following 
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facts taken from Scofield and Oliver and Woodley et-al.^^^^ regarding rain 
clouds : 

1. Bright clouds in the visible imagery produce more rain than 
darker clouds. 

2. Bright clouds in the visible imagery and clouds with cold tops 
in the infrared imagery which are expanding in areal coverage 
produce more rainfall than those not expanding. 

3. Decaying clouds produ,ce little or no rainfall. 

4. Clouds with cold tops in the IR imagery produce more rain than 
those with warmer tops. 

5. Clouds with cold tops that are becoming warmer produce little or 
no rain. 

6. Merging of cianulonimbus clouds increases the rainfall rate of the 
merging clouds. 

7. Most of the significant rainfall occurs in the upwind Cat anvil 
level) portion of a convective system. The highest and coldest 
clouds form where the thunderstorms are most vigorous and the rain 
heaviest. These cold clouds get thinner downwind and look warmer 
in IR imagery as the anvil material blows away from its origin 
over the updraft. 

The Scofield/Oliver technique was tentatively selected as a model for making 
daily estimates of precipitation over Virginia because of its great success 
in estimating flood-producing rains over the United States in near real time 
Scofield and Oliver^^^^, Scofield . Experiment with this technique, 
which was developed primarily for very heavy convective rains during the warmer 
months of the year, suggested that modifications would be advisable when 
applying it to light rain situations, to synoptic scale storms, and 
(especially) to winter storms. Time being of the essence in the operational 
mode, it was decided to follow the hourly changes and movements of storm clouds 
in the enhanced infrared (EIR) only, rather than the half-hourly in both the 
EIR and visible imagery, as called for in the Scofield/Oliver scheme. 

Scofield tried this approach on the large synoptic storm over the north- 
eastern United States on October 9, 1976, with encouraging results. This 
variation of the Scofield/Oliver technique will be referred to hereafter as 
the S technique. 
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Hourly radar charts from the National Facsimile Circuit of NWS, and hourly 
radar reports from Patuxent, Maryland and Bristol, Tennessee were used to 
locate areas of heavier rain, but (generally speaking) not to estimate amounts. 
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SECTION 2 . PRELIMINARY STUDIES 


Seven storms occurring between May and November 1978 were analyzed by 
hand. Isohyets of hourly rainfall estimates, based on the S technique 
slightly modified by radar reports, were drawn on the map of Virginia and 
its immediate vicinity. From these isohyetal charts, estimates for 2S 
stations equipped with rain gages were tabulated for each hour from 1200Z 
Cnoon Greenwich Mean Time) on one day to 1200Z the following day. The hourly 
estimates were summed and compared with observed rainfall amounts at the 2S 
stations for the same time period. The results ranged from excellent to poor. 
Obviously, to be operationally effective, estimates adjacent to these stations 
would have to be adjusted to the gage readings. Furthermore, the 25 stations, 
which report daily on National Facsimile and/or Service C Teletype Circuits, 
should be supplemented by as many volunteer stations as feasible. 

Eventually a simplified and rapid method of estimating rain from satellite 
imagery was designed as a possible alternative to the more laborious and 
time-consuming S technique. This quick (Q) technique assigned rain for a 
given hour in accordance with the cloud top temperature in the EIR imagery 
at the end of that hour. (Later the cloud top temperature at the beginning 
of the hour was used.) The digital enhancement (MIb) curve was used: areas 
shaded medium gray (-32 to -41 ®C) were assigned 0.05 inch of rain; areas 
shaded light gray (-41 to -SZ^C) were assigned 0.10 inch, areas of dark gray 
C-52 to -58®C) were assigned 0.25 inchj black, areas (-58 to -62“C) were 
designated 0.50 inch; repeat gray level areas (-62 to -80®C) were designated 
1.00 inch, and white areas (below -80®C) were assigned 1.50 inches. (See 
Table 2-1.) 


Table 2-1. Cloud Top Temperature vs Precipitation Used in 1978 Storms 


Enhanced Shade 
(Mb curve) 

Medium gray 

Light gray 

Dark gray 

Black 

Repeat gray 
White 


Temperature Range 
C“C) 

-32 to -41 
-41 to -52 
-52 to -58 
-58 to -62 
-62 to -80 
Below -80 


Precipitation 

(inches) 

0.05 

0.10 

0.25 

0.50 

1,00 

1.50 
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Th« Q technique gave larger errors than the S technique in estinuting precipi^ 
tation at the 2S daily reporting stations for the seven 1973 storms. However, 
it was never intended that either method be used alone as an operational tool, 
Ideally, any operational satellite technique should be corx'ected for "ground 
truth" observed rain at all gages in the network, and modified by interpolation 
of the gage readings. Or, put another way, an isohyetal field based on the 
gage readings alone should be modified by reference to the satellite estimate 
field. Satellite methods have the advantage of covering the entire field, but 
fail to provide exactly accurate estimates at given points. On the other hand, 
gage networks provide exact measurements at fixed points, but the interpolations 
are inexact. This suggests that optimal results are likely to be obtained by 
a judicious blend of the two. With this in mind, nine methods of estimating 
24-hour precipitation over the state were designed and tested on the seven 
storms of 1973. These nine methods will be referred to as G, S, Q, Sm, Qm, Sa, 
Qa, Sam and Qam. 


Method G used 24-hour precipitation amounts from a network of 63 rain gages 
to interpolate the estimate field. The interpolation method is described in 
the next section. This network consisted of 25 stations which report daily 
on the NWS facsimile and teletype circtuts, and 38 stations reporting several 
months after the fact in NWS Climatological Data Summaries CCDS) . Forty-nine 
stations’ reports from the CDS were held in reserve as control data for 
verification purposes. CSee Figure 2.1) 


The S and Q methods are the S and Q techniques described above. The Sm method 
produces an estimate field from the mean of the S value and G value at each 
point in the field. The Qm method takes the mean between the Q and G values 
at each point to obtain the Qm field. 


Considering the gage reading to be ground truth at the gage location, the 
error for S at each analysis gage location was computed and an error field 
derived by interpolation. S values at all points in the state were adjusted 
by this error field to get the Sa field. This is the Sa method. The same 
procedure was used to adjust the Q values and obtain the Qa field. This is 
the Qa method. 

The Sam method takes the mean between the Sa and G values at each point to 
obtain the Stsm field. The Qam method obtains the Qam field from the mean 
of the Qa value and G value at each point. (In actual practice the inter- 
polations were made only for the control station locations, in order to verify 
each method.) 
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Each o£ tha nina bmcHs!49 wu tastadi for each of tha Sevan storms in 1978, 
using tha 49 control stations for verification. Table 2-2 gives composite 
scores for each mathod for the seven storms taken together, including mean 
absolute error, mean algebraic error (bias), and tha correlation coefficient 
for each method versus obsairved precipitation. The Qa method obtained the 
bast scores, but the Qam and Sara, ranking second and third, respectively, 
ware almost as good. Tha 9 method ranked a poor last, and tha S method Is 
next to the last, but decidedly superior to Q. The G method is in the middle 
of tha group, but it shows tha highest correlation coefficient on three of 
tha storms, and the least absolute error on two stoiws. Table 2-3 through 
2-9 show the performance of each method on each individual storm. 

These scores suggest that although either satellite method taken alone will 
seldom produce the best estimates state-wide, estimates based on the Q method, 
when adjusted by observed gage readings, may frequently be more accurate than 
estimates based on interpolations between observed gage measurements alone. 

The scores also indicate that it is safe to adopt the various Q methods while 
abandoning the labor-intensive, time-consuming S methods. 
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Tabl* 2-2* Compositt Scor«» fpy each Method for Seven 1978 Storms Combined 


Method 

Correlation 

Coefficient 

Mean Error 
(Absolute) ♦ 

Mean Error 
(Algebraic)* 

G 

0.824 

(63 

0.200 

(4) 

-0.014 (3) 

S 

0.713 

(3) 

0,257 

(8) 

-0.123 (8) 

Q 

0.481 

(9) 

0.336 

(9) 

-0.197 (9) 

Sm 

0.S3S 

(4) 

0,200 

(4) 

o 

o 

o 

1 

Qm 

0.820 

(73 

0,235 

(7) 

-0 109 (7) 

Sa 

0.827 

(5) 

0.209 

(6) 

-0.016 (5) 

Qa 

0.84S 

(13 

0.194 

(1) 

-0,006 (1) 

Sam 

0.341 

(23 

0.196 

(3) 

-0,014 (3) 

Qam 

0.841 

(23 

0.194 

Cl) 

-0,010 (2) 

N»317. 

Mean Rain at 

Control 

Stations : 

0.574 

inch. 

NOTE; 

Numbers in 

parentheses indicate rank. 

♦Inches. 


Table 

2-3. Scores for Storm May 4-5, 1973 

Method 

Correlation 

Coefficient 

Mean Error 
(Absolute)* 

Mean Error 
(Algebraic) * 

G 

0.915 


0.188 


-0.01 

S 

0,598 


0.441 


-0.36 

Q 

-0.192 


0.773 


-0.69 

Sm 

0.875 


0.264 


-0.19 

Qm 

0.835 


0.420 


-0.38 

Sa 

0.399 


0.213 


-0.06 

Qa 

0.904 


0.198 


0.00 

Sam 

0.914 


0.194 


-0.04 

Qam 

0,913 


0.187 


-0.01 

N«46. 

Mean rain at 

control 

stations : 

1.370 inches. 


NOTE: Best scores are underlined. * Inches. 



TabI«2-4. Scores for Storm May 31-June 1, 1978 


Method 

Correlation 

Mean Error 

Mean Error 


Coefficient 

(Absolute)* 

(Algebraic) 

G 

0.776 

0.0S2 

+0.04 

S 

0.838 

0.015 

+0.01 

Q 

0.670 

0.015 

0.00 

Sn 

0.907 

0.023 

+0.03 

Qm 

0.828 

0.026 

+0.02 

Sa 

0.73S 

0.044 

+0.03 

Qa 

0.314 

0.050 

+0.04 

Sam 

0.767 

0.048 

+C.04 

Qam 

O.SOS 

0.051 

+0.04 

N-4S. 

Mean rain at control 

stations; 0.014 

inch. 

NOTE; 

Best scores underlined. * Inches. 



Table 2<^S. Scores for Storm June 

8-9, 1978 

Method 

Correlation 

Mean Error 

Mean Error 


Coefficient 

(Absolute)* 

(Algebraic) 

G 

0.867 

0.130 

-0.03 

S 

0,550 

0.259 

-0 . 10 

Q 

0.142 

0.285 

0.00 

Sm 

0.826 

0.149 

-0.06 

Qm 

0.754 

0.170 

-0.01 

Sa 

0.313 

0.171 

*0 . 03 

Qa 

0.810 

0.150 

-0.03 

Sam 

0.862 

0.139 

0.00 

Qam 

0.849 

0.132 

-0.03 

N«44. 

Mean rain at control 

stations; 0.459 

inch. 


NOTE; Best scores underlined. * Inches. 
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Table 2-6, Scores for Stotm June 21-22, 1971^^ 


Method 

Correlation 

Coefficient 

Mean Error 
(Absolute ) * 

Mean Error 
(Algebraic) 

G 

0.363 

0.301 

♦0.04 

S 

0.282 

0,314 

-0.03 

Q 

0.43S 

0.244 

-0.06 

Sm 

0.368 

0.278 

0.00 

Qm 

0.439 

0.263 

-0.01 

Sa 

0.409 

0.320 

♦0.02 

Qa 

0.S41 

0.2S4 

♦0.04 

Sam 

0.406 

0.295 

♦0 , 03 

Qam 

0.477 

0.276 

♦0.04 

N»46. 

Mean rain at control 

stations; 0.369 

inch. 


NOTE; Best scoros underlined. * Inches. 


Table 2-7. Scores for Storm July 14-lS, 197S 


Method 

Correlation 

Coefficient 

Mean Error 
(Absolute)* 

Mean Error 
(Algebraic) 

G 

0.576 

0.192 

+0.05 

S 

0.42S 

0.138 

-0.02 

Q 

0.280 

0.193 

-0.08 

Sm 

0.595 

0.174 

♦O.OI 

Qm 

0.581 

0.165 

-0.02 

Sa 

0.599 

0.173 

+0.03 

Qa 

0,544 

0.190 

♦0.08 

Sam 

0.632 

0.176 

♦0.04 

Qam 

0.567 

0.188 

♦0.07 

N»43. 

Mean rate at control 

stations: 0.222 

inch. 


NOTE: Best scores underlined. *Inches. 


2-7 



Tablt 2-8, Scorts for S^ozm July 31-August I, 1973 


Method 

Correlation 

Coefficient 

Mean Error 
(Absolute) * 

Mean Error 
(Algebraic)* 

Q 

0.19S 

0.415 

-0.20 

S 

0.499 

0.384 

-0.26 

Q 

0.366 

0.403 

-0.15 

Sn 

0.522 

0.369 

-0.23 

Qn 

0,462 

0,365 

-0.17 

Sa 

0.438 

0.380 

-0,18 

Qa 

0,453 

0.334 

-0.17 

San 

0.411 

0.338 

-0.19 

Qan 

0.3S7 

0.395 

-0.19 

N-48. 

NOTE: 

Mean rate at control stations: 0.S36 inch. 

Best scores underlined. * Inches. 

■ 

Table 2-9, Scores for Storm November 16-17, 1978 

Method 

Correlation 

Coefficient 

.Mean Error 
(Absolute)* 

Mean Error 
(Algebraic)* 

G 

0.351 

0.105 

+0.01 

S 

0.109 

0,180 

-0.10 

Q 

-0.440 

0.418 

-0.40 

Sm 

0.638 

O.llS 

-0.05 

Qm 

0.657 

0.220 

-0.19 

Sa 

0,574 

0.144 

+0,02 

Qa 

0.302 

0.123 

0.00 

San 

0.7SS 

0.114 

+0.02 

Qan 

0.846 

O.llO 

+0.01 


N"45. Mean rate at control stations; 0.SS4 inch. 
NOTEi Best scores underlined. ^Inches. 
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SECTION 3 ~ DESCRIPTION OF AUTOMATED TECHNIQUE 

Photographs of Virginia and vicinity taken by the geostationary satellite 
located' over the equator at 75* west longitude consist of very small picture 
elements (pixels) . Due to foreshortening at the latitude of Virginia, each 
pixel is a photograph of a rectangular area 9.7 kilometers north=^south and 
6.2 kilometers east-west. There are 4879 of these cells in the rectangle en- 
closing the state of Virginia and its immediate vicinity. The information 
in each pixel is the cell's cloud top temperature to the nearest whole degree 
Celsius. If digital cloud data could be accessed at need, a linear relationship 
between cloud top temperature and hourly precipitation might be used. This 
would be a great refinement of the scale shown in Table 2-1, 

NASA was fortunate in arranging for receipt of the digital IR and visible 
imagery for Virginia on a continuing basis. The dedicate source is the 
National Meteorological Center (NMC) of NWS located at Camp Springs, Maryland. 
Hourly digital imagery, both IR and visible, is stored in NMC's data tank. 

Each picture, if needed, must be acquired by NASA (later by VPI5SU) via auto- 
dial within 18 hours of picture-taking time (before it is erased in order to 
store more recent pictures.) 

As a consequence of this data breakthrough, a number of options and strategies 
under consideration were abandoned. These included time-consuming hand 
analysis of satellite imagery by a professional meteorologist, transfer of 
hourly isohyetal analysis to bit-pad equipment using a cursor, and complicated 
summation of the isohyetal fields to get 24-hour totals for all cells. Since 
these slower procedures were abandoned, no description is necessary. 

The system now in use consists of two computer programs written in FORTRAN. 

See Appendix A for a complete listing of these programs. The first program 
computes the satellite estimate (Q) field from the hourly digital IR imagery, 

The second program computes the final estimate for the entire state area by 
comparing five preliminary estimates of 24-hour precipitation (Q, G, Qa, Qm 
and Qam) with control raingage readings and determining which of the five 
methods gives the best estimate for the day. The final estimate is then 
produced by incorporating control gage readings into the winning method. 

The satellite estimation program consists of the main drive routine RAINDR 
and five subroutines: PREPRO, RDATA, QUICK, LOOKUP, and RAINSM (see Appendix A). 

ORlGimi 
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l§ 
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RAINDR processes each hour of satellite IR imagery as follows: First, in sub- 

routine PREPRO, the raw data (as received from NMC via autodial) is stripped 
down to the minimum sire which still covers the entire state. Next, each 
pixel in the reduced data set is converted from temperature to rainfall amount 
in subroutine QUICK (the Q method) using the temperature to rainfall conversion 
table generated in subroutine LOOKUP. This new data set is then added to the 
total rainfall data set in subroutine RAINSM. After all 24 hours of IR imagery 
have been processed the sum in the total rainfall file is the 24-hour satellite 
estimate - the Q method estimate for the day. 

The final estimation program consists of the main driver routine MAIN and 7 
subroutines: GNDTRU, FILSIA, SORT, INTERP, NXTPNT, RANGE, and WEIGHT. The 

final rainfall estimate is computed by MAIN as follows: The satellite estimate 

computed previously in the first program, and the gage readings for each 
reporting ground station are entered in subroutine FILSTA. Subroutine INTERP 
then computes the interpolated rainfall and error values for each control 
station in the system using these values at each analysis station. Subroutines 
SORT, NXTPNT, RANGE, and IVEIGHT are used in the interpolation scheme which 
will be described later. Once the interpolated rain and. error values have 
been computed for each control station, the five preliminary estimates are 
computed. G is simply the interpolated rainfall field; Q is the satellite 
estimate computed previously; Qa is the satellite estimate (Q) plus the inter- 
polated error value (the error value adjusts the satellite estimate); Qm is 
the mean of the interpolated rainfall field and the satellite estimate; and 
Qam is the mean of Qa and the interpolated rainfall field. From these 
values an error is computed for each of the five estimates and the absolute 
value summed for all control stations. The method with the lowest total error 
is proclaimed the official method or "winner" for the day. Once the winner 
has been found, the control station gage readings are incorporated into the 
winning method and the final estimate is computed. When G wins, subroutine 
INTERP computes the interpolated rainfall field for the entire state using 
both control and analysis stations. When Q wins, the previously computed 
satellite estimate is the field estimate (but all observed gage readings 
override the satellite in their cells.) IVhen Qa wins, INTERP computes the 
interpolated error field for the entire state, which is then added to the 
satellite estimate, producing the final estimate. When Qm wins, INTERP 
computes the interpolated rainfall field and the final estimate is the mean 
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of this field and the satellite estimate. When Qam wins, INTERP computes 
both the interpolated rain and error fields, the satellite estimate is added 
to the' error field (Qa field) and the final estimate is the mean of this 
field and the interpolated rain field. See Figure 3-1 for general flow diagrams. 

The curve first used in automated technique relating cloup top temperature to 
hourly precipitation actually consists of two linear curves with different 
slopes. -32“C corresponds to 0.01 inch and -56*C corresponds to 0.25 inch; 
that is, rain increases by 0.01 inch for each l“C decrease in temperature. 
Thereafter rain increases by O.OS inch for each 1“C decrease in temperature 
of the cloud top. Thus -S3“C corresponds to 1.60 inches of rain per hour, 
and the atmosphere should seldom become colder than this. 

For the storm of August 15-16, 1980, the first storm analyzed by the new 
automated technique, a number of stations which report hourly precipitation 
amounts were compared with the corresponding cloud top temperature . in the 
pictures at the beginning and end of the hour. A new curve was generated 
from the resulting scatter diagram. For this curve, -24®C corresponds to 
0.01 inch, and -56®C corresponds to 0,17 inch, for an increase of 0.01 inch 
for each 2“C decrease in temperature. Thereafter rain increases 0.06 inch 
for each 1®C decrease in temperature, reaching 1.80 inches at -83®C. (See 
Figure 3-2. ) 

Both curves have been applied to each of the 1980-81 storms studied. The 
verifications are listed in Tables 5-1 through 5-22. As a result, the old 
curve has been adopted as officials 

The stations used in analysis and verification of 1978 storms (Figure 2-L) 
are not uniformly distributed, but these were the only available stations 
which report at or near 7 a.m. EST (1200Z.) By August 1980 several of these 
stations had been discontinued, while a considerable number had been added. 

Figure 3-3 shows the location of stations used in analyzing and verifying 
the storms of 1980 and 1981. Twelve Nationwide Agricultural Touchtone 
System C^ATS) stations, recruited by VPI5SU, are included. The distribution 
still leaves eastern Virginia and the southwestern counties under populated. 

These gaps are to be filled with NATS stations in the near future. Also the 
delayed stations, which report in the CDS several months late, will have to 
be replaced unless the observer can be induced to report in real time via NATS. 
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Figure 3-1, Program MAIN 
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Figure 3-1. Lowest Total Error 
3-6 
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Th« automated technique generates five fields of 24-hour precipitation 
estiaates for the state of Virginia « tests each field against control data, 
and chooses (,foT the da/ in question) the field which makes the best score. 
The fields I designated G, Q, Qm, Qa and Qan, ma/ be subdivided into prelimi- 
nary and final fields as the technique proceeds. Each field consists of 2885 
rectangular cells which cover the state. 


All 24-hour rain gage readings taken at or near 7 a.m. EST (noon GMT) are 
ingested by the mini -computer, which interpolates between the analysis 
stations to get the preliminary G field. The control gage readings are stored 
by the computer for future use. The preliminary G field provides a 24-hour 
estimate of precipitation for each cell in the state for the period from 7 a.m. 
to 7 a.m. EST C1200Z to 1200Z.) 

Hourly digital IR imagery for the 24-hour period is received via auto-dial 
from the NMC data tank. The mini-computer assigns an hourly precipitation 
estimate to each cell in Virginia, using the curve relating cloud top 
temperature fig precipitation. The hourly precipitation estimates are summed 
by the coraptuer to provide a 24-hour total for each cell. This is the pre- 
liminary Q field. 


Each observed gage measurement is considered ground truth for the 9.7 by 6.2 
km cell in which the gage is located. This observed value overrides any 
estimate based on the satellite data in that cell. Therefore the computer 
compares the observed amount at each analysis station with the satellite CQ) 
estimate for that cell, the difference being the error in the Q field in 
that cell. An error field for Q is then obtained by interpolation between 
all analysis stations, but not the control stations. This error field is 
used to modify all cells in the Q field. These modified cells comprise the 
Qa field, or the adjusted satellite field. 

FPr example, suppose that cell A lies somewhere between cells B, C, D and E, 
which contain analysis stations (see Figure 3-4.) To obtain the probable 
error in the Q field at cell A, the Q field errors at B, C, D and E are 
given weights in proportion to the inverse (reciprocal) of their respective 
distances from A. Distances between cells are measured from cell center to 
cell center. Say that B lies SO km southwest of A; C lies 40 km westnorth- 
west of A; D is 70 km north of A; and E is 60 km east of A. We'll asr-ime 
that no other analysis station lies within 100 km of A (the technique's 
arbitrary cut-off distance for influence.) The influence B exerts on A 
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Q ■ 1.40 
O0G « 0.30 
Qa^ 0.80 



B 0G « 1.80 


Q « 0.88 
E-0 Q « 0.88 
Qa ■ 0.88 


Qa - 1.80 

Figure 3>S. Interpolation Procedure as used to Derive G Field 


is proportional to the inverse of the distance from A to B; that is, the re- 
ciprocal of 50, or 0.02. C's influence is the reciprocal of 40, or 0.02S. 

D's influence Is 0.014, and E's influence is 0.017. The total influence is 
their sum, wh^ch is 0.076. This assigns B 26% of the weight, C 33%, D 19% 
and E 22%. If the Q error at B is, say, -0.97, at C is -0.30 inch, at D is 
4-0.60 inch, and at E is zero, then B will contribute -0.2S2 inch to the error 
at A, C will contribute -0.099 inch, 0 will contribute 4-0.114 inch, and E 
will contribute zero. The algebraic sum is -0.237, or -0.24 inch, the pro- 
bable error in the Q field at A. If the Q value at A is, say, 1.12 inches, 
then the Qa value at A will be 1.36 inches. 

At this point the computer could perform this adjustment of the Q estimates 
for each cell in the state to derive the Qa field, but in practice it now 
makes the adjustment at the control station only. Note that the Qa value at 
each analysis station will equal the observed gage measurement, but that Qa 
will rarely equal the gage measurement at each control station. 

The next step is to derive the Qm field. The computer performs this task by 
taking the mean of G and Q values in each cell. However, at each analysis 
station the gage reading will override any differing Q value in that cell. 
That is, in each analysis gage cell, G*Q*>Qm. 

The Qam field is derived in the same way, by determining the mean between the 
G and Qa values in each cell in that state. 

For verification, an error is computed for each field at each of the control 
stations. The errors for each field are summed, and the field with the least 
absolute error is chosen as the "best” precipitation field for the 24-hour 
period in question. This is sti a preliminary field. It must now be 
adjusted by means of a control gage readings, which up to this point have 
not been used in any of the analyses, to get the final official precipitation 
field for the day. 


The example in Figure 3-4 is further expounded in Figure 5-5. Hypothetical G 
and Q values have been entered at points B, C, D and E, and are compatible 
with the errors in Q shown in Figure 3-4. Multiplying the G values, i.e. 

1,60, 1.40, 0.80 and 0.88 by 26%, 33%, 19% and 22%, respectively, and 
summing, we get an interpolated G estimate, G^^, of 1.22 inches at A. The 
Qm estimate at A is the mean of the G and Q values, or 1,17 inches. The Qam 
estimate is the mean of the G and Qa values, or 1.29 inches. (Note that G 
and Qa do not have the same value at A, since G is merely interpolated. 
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not ground truth. Also note that Q at ‘point A is not an interpolation of 
the Q values at B, C, D and E, since it has been derived directly from, the 
satellite imagery.) 

Assume that point A is a control station. In this event, the measured amount 
at A is ground truth and will be used to determine the error at A of each 
estimate field. The estimates, of course, will have no influence on the 
observed amount. In our example, if the observed precipitation at A is 1.32 
inches, then the errors areG«-0.10, Q*-0.20, Qm»-0.,1S, Qa«+0.04, and Qam*-0.03 
inch. 

The computer compares the estimate obtained by each of the five methods at 
each control point with the corresponding gage measurement and computes the 
correlation coefficient, R, for each method. Normally the method with the 
smallest total absolute error will have the highest correlation coefficient, 
but when two methods are virtually equal in least total absolute error, R 
may be used as a tie-breaker. 


■ POOR QUALlfY 
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SECTION 4 - CRITIQUE OF AUTOMATED TECHNIQUE 


The temperature of cloud tops depicted in infrared satellite imagery accounts 
for some of the variability in precipitation under those clouds. A scheme 
which attempts to estimate precipitation from cloud top temperature alone 
suffers from great over-simplification and can attain only limited success. 
Alternatively, the inclusion of other parameters (not including observed rain- 
fall) complicates, the procedure, necessitating either more complex computer 
modelling or the services of an on-the-job satellite meteorologist. The 
meteorologist would take cognizance of the following rules, among others; 

Cl) Shields of cirrus and cirrostratus are cold, but taken alone 
produce no precipitation. 

(.2) In the convective regimes of warm weather, thunderstorms and heavy 
rain tend to be concentrated in the sharply defined upwind (at cirrus level) 
portion of the cold cloud mass rather than in the filmy ill-defined downwind 
portion, 

C3) Convective cells are often embedded in synoptic-scale systems, 
particularly in colder months, but are detected only with difficulty by the 
satellite meteorologist. 

(.4) Heavier precipitation generally occurs in the comma cloud, in areas 
of maximum positive vorticity advection, and near the inflection point of the 
jet stream over or near a front as delineated in the satellite imagery. 

C5) On occasion the coldest clouds over the state will be non-precipi- 
tating debris from old systems, while somewhat warmer clouds are producing 
copious rain or snow. 

The expertise of a trained satellite meteorologist on a continuing basis is 
required to recognize and cope with situations such as those referred to in 
the above examples. It is with this in mind that we should evaluate any 
shortcomings of the automated technique, which does have the decided advan- 
tage of requiring no expert input. 

Missing data, and garbled data received by auto-dial from the NMC computer, 
have plagued the investigation and may well plague operations in future. 
CFifteen out of 24 hours of data were missing during the storm of March. 5-6, 
1981, and 14 hours were missed during the September 17-18, 1980 storm. 

(see Table 4-1) Interpolation algorithms are fairly reliable when a single 
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Table 4>li : Pixel Shift, Observed Storm Total C^recipitation) , and Number 

of Control Stations, Analysis Stations and Missing Hours of 
Satellite Data for Storms of 1980 and 1981 


Storm Date 

Pixel 

Shift 

Storm 

Total# 

Control 

Stations 

Analysis 

Stations 

Missing 

Hours 

1980 

Aug lS-16 

3S IE 

7.61 

47 

63 

11 

Sep 9-10 

3S 2E 

4.75 

SI 

61 

12 

Sep 16-17 

2S 0 

0.24 

52 

66 

10 

Sep 17-18 

SS IE 

20.36 

56 

73 

14 

Sep 24-25 

2S 2W 

27.12 

56 

69 

10 

Sep 25-26 

2S IN 

3,61 

S3 

68 

7 

Oct 24-25 

0 IE 

41.13 

49 

63 

3 

Nov 15-16 

4S 0 

5,71 

48 

59 

5 

Nov 17-18 

IS IW 

39.75 

52 

71 

4 

Dec 9-10 

2S 0* 

11.24 

53 

80 

7 

1981 

Mar 4-5 

IS 0 

25.56 

57 

6S 

6 

Mar 5-6 

1S2W 

2.23 

54 

68 

15 

Mar 16-17 

2S 0* 

4.85 

50 

62 

7 

Mar 22-23 

2S 0* 

25.26 

S3 

67 

12 

Mar 29-30 

2S 0* 

9.99 

51 

64 

4 

Mar 30-31 

IS 2W 

22.20 

51 

58 

7 

Apr 14-15 

2S 0* 

12.40 

53 

74 

6 

May 15-16 

2S 0* 

5.08 

53 

62 

7 

May 18-19 

2S 0* 

31.66 

53 

66 

9 

May 19-20 

2S 0* 

26.47 

53 

64 

9 

Jul 21-22 

IN 2W 

2.69 

17 

32 

6 


# Inches. *Overall mean shift. 
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hour of data is missing* A particularly unsatisfactory situation is sdvnral 
hours of missing data at the very beginning or end of the day, with no anchor 
point for the interpolation, Recently reception has improved, but three ox' 
four missing hours per day is not uncommon. Garbling, a very serious problem 
in the early stages of the investigation, has been nearly ollminated. 

The digital imagery received via auto-dial consists Of an array of 129 by 129 
pixels, the picture center (the 65th pixel of the 6Sth row) is supposed to 
be at 3?" north latitude, 79** west longitude, about five kilometers southwest 
of Brookneal, Virginia. Cloud features in actual satellite BIR photographs 
(Figure 4-1) wore carefully compared with the corresponding features in the 
digital Imagery accessed by auto-dial (Figure 4-2.) In 13 storms occurring 
between August 1980 and March 1981, geographical location of the Clouds was 
somewhat different in the two typos of imagery. In general It was necessary 
to shift the digital pixels to the north to bring cloud patterns into line 
with their counter{)arts in the EIR photographs. Put another way, it was 
necessary to examine pixels to the south of the 65th pixel of the 65th row 
to find clouds over 37° noi'th, 79° west (picture center). The actual shift 
of the stations varied from storm to storm, and even from hour to hour, but 
averaged two pixels south and almost one-half pixel east for all cases. 
n»at is, instead of looking for the cloup top temperature over Richmond at 
the 88th pixel in the 59th row (the nominal position for Richmond) we should 
expect it in the 88th pixel (or perhaps the 89th) in the 6lst x*ow. We have 
not found the reason for this shift, but NMC is aware of it and is attempting 
to solve the problem. 

Figure 4-1 is the BIR photograph taken by the geostationary satellite at 
0300Z August 16, 1980, showing cloud features in relation to geographical 
landmarks. Figure 4-2 is the corresponding digital imagery for 0300Z August 
16, 1980, accessed by auto-dial. Careful examination of clouds in relation 
to landmarks shows that corrections in registration must bo made in the 
digital product to bring it into alignment with the EIR photograph (which has 
been 'quite accurately rogistoired by the National Environmental Satellite 
Service). These corrections are indicated in Table 4-2. T1\e shift column 
shows the direction of the error. The correction must bo made in the opposite 
direction. 
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Figure 4-1. EIR Photograph taken at 0300Z August 16, 1980 

showing Cloud Features in Relation to Landmarks 
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Tabl« 4-2. Apparent Shift in Landmark Positions in Digital 
'laagery for OSOOZ August 16, 1980 


Landmark 

Nominal 

Row 

Position 

Column 

Observed Position 
kow Column 

Shift 

Southeast comer 
of Virginia 

70 

109 

73 

108.5 

3S 0.5W 

Cape Charles 

63 

108 

65. S‘ 

107 

2.5S IW 

Corner west of 
Charlotte, NC 

85 

35 

87 

38 

2S 3E 

Comer west of 
Laurenburg, NC 

90 

55 

92.5 

56 

2.SS IE 

Southwest comer 
of Delaware 

48 

112 

48 

no 

0 2W 

Mean Shift 





2S O.IE 

Similar shifts were determined 

for the imagery 

for 2000Z, 21002 and 

22002 of 

August IS, 1980, 

and for OOOOZ, 

OlOOZ, 0200Z, 

04002 and 

0600Z of August 16, 1980 


The average shift for the storm day was 3 south 1 east. Four additional land- 
marks were used in the full day analysis: the southwest tip of Maryland (row 

39/column 58); the junction of Ohio, Kentucky and West Virginia (row 48/colunm 
13); the junction of Virginia, Tennessee and North Carolina (row 69/column 26); 
and Picture Center (row 65/colunm 65) . 

Another source of possible error, which applies equally to the EIR photographs 
and the digital imagery, is parallax due to the angle at whidi the satellite 
over the equator at 75^ west longitude photographs the clouds, A 40,000 foot 
cloud top over Virginia appears to be about 6.5 nautical miles north of its 
true position. A 30,000 foot cloud top has an apparent displacement of 5.5 
nautical miles, or 9.8 km, the north-south dimension of a pi.xel in the latitude 
of Virginia. The computer program carries a correction for this parallax 
phenomenon. 


SECTION 5 - ANALYSIS OP EXPERIMENTAL STORMS OF 1930-81 


The automated technique ^as tested on 21 storms which occurred in 1980 and 
1981. Composite scores for these storms are shown in Table S-1. The scores 
attained by each of the five methods are tabulated in Tables s-2 through 5-22. 
Table 5-23 shows how frequently each method had the highest correlation co- 
efficient (estimated vs observed precipitation) , the least mean absolute error, 
and the least mean algebraic error (overestimates and underestimates). The G 
method makes the highest scores overall, somewhat better with the old curve 
(cloud top temperature vs precipitation) than with the new curve. The Qara 
method does well in least mean absolute error, particularly with the new curve. 
Qa scores well in bias (mean algebraic error), leading on six storms (29 
percent) with each curve. 

IVhile the technique is designed primarily to estimate growing season rain, 
most of it convective in nature, only five of the 21 storms in the study were 
typical airmass storms. The remaining 16 were associated with fronts and 
synoptic-scale low pressure systems. Appendix B gives a brief description of 
the synoptic situation for each of the storms. 

The original computer program was completed and the link with NMC’s computer 
was established in August 1980. Four summer storms were accessed and analyzed 
between August 15 and September 21, 1980. Among these the ’’storm" of 
September 16-17 produced very little rain. Subsequently six storms in autumn, 
three in winter and seven in spring were studied. Storms over Virginia were 
scarce in the summer of 1981, and the storm of July 21-22, 1981 alone was 
accessed and analyzed. Only 17 control stations and 29 analysis stations were 
available for this storm; the results therefore should be viewed with caution. 

In view of the scarcity of summer airmass storms during the experimental 
period from .August 1980 to July 1981, the uniform rainfall patterns observed 
should not be surprising. A good example of uniform distribution is the 
weak storm of April 14-15, 1981. (See Figure 5-1..) Heaviest rain in Virginia 
(one-half to three-quarters of an inch) occurred in an eliptical area 
extending between Trout Dale and Hillsville in the southwest part of the 
state. Most other sections had between 0.1 and 0.4 inch, and only three 
reporting points had no rain. I'Tneneyer such uniformity prevails the inter- 
polation algorithm incorporated into the program makes the "gages only" 
routine hard to beat by any satellite method. 
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Figure 5-1. Precipitation Pattern on April 14-15, 1981 
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Table 5-1 


Conposite Scores for each Method for 21 Storms in 1980-81 


Old Curve 

Q 

G 

Qa 

Qm 

Qam 

Mean Absolute 
Error (inches) 

0.525 

0.123 

0.181 

0.295 

0.145 

Mean Algebraic 
Error (inches) 

*0.310 

+0.004 

+0.043 

+0.157 

+0.024 

New Curve 






Mean Absolute 
Error Cinches) 

0,547 

0,123 

0.172 

0.311 

0.143 

Mean Algebraic 
Error (inches) 

+0.428 

+0.004 

+0,042 

+0.216 

+0.023 


N»1062. Best scores underlined. 
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Table S-2. 

Scores for Storm August IS- 

16, 1980 


0l<4 Curve 

Q 

G 

Q» 

Qb 

Qan 

Correlation 

Coefficient 

0.687 

0.836 

0.800 

0.768 

0.863 

Keen Absolute 
Error* 

0.2S9 

0.093 

0.110 

0.1S9 

0.089 

Mean Algebraic 
Error* 

* 0.206 

-0.001 

+0.037 

+0.103 

+0.013 

New Curve 






Correlation 

Coefficient 

0.71S 

0.S36 

0,771 

. 0.788 

0.853 

Mean Absolute 
Error* 

0.411 

0.093 

0.116 

0.229 

0.091 

Mean Algebraic 
Error* 

+0.403 

-0,002 

+0.02$ 

+0 . 201 

+0.010 

N«47. Mean 

Rain at Control Station: 

0,162 inch. 

* Inches . 


Table S-3. 

Scores for Storm September 9 

-10, 1980 


Old Curve 

Q 

G 

Qa 

Qm 

Qam 

Correlation 

Coefficient 

0.466 

0.832 

0.624 

0.687 

0.770 

Mean Absolute 
Error* 

0.214 

0.056 

0.089 

0.12S 

0,065 

Mean Algebraic 
Error* 

+0.170 

-0.020 

+0.001 

+0.075 

-0.010 


New Curve 


Correlation 

Coefficient 

0.490 

0.832 

0.647 

0.718 

0.773 

Mean Absolute 
Error* 

0.240 

0.056 

0.077 

0.141 

0.060 

Mean Algebraic 
Error* 

+0.201 

-0.020 

-0.007 

+0.090 

-0.014 

N-Sl. Mean 

Rain at Control 

Stations : 

0,093 inch. 

* Inches. 
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Tabl« 5«4. Scores for Stom September 16-17, 1980 


Old Curve 

Q 

G 

Q* 

Qm 

Qam 

Correlation 

Coefficient 

-0.087 

0.4S2 

0.077 

-0.020 

0.193 

Mean Absolute 
Error* 

0.024 

0.006 

0.013 

0.015 

0.010 

Mean Algebraic 
Error* 

+0.01S 

-0.002 

+0.004 

+0.006 

+0.001 

New Curve 






Correlation 

Coefficient 

-0.161 

0.452 

-0.086 

-0.118 

0.075 

Mean Absolute 
Error* 

0.072 

0.006 

0.014 

0.039 

0.010 

Mean Algebraic 
Error* 

+0.065 

-0.002 

+0.005 

+0.032 

+0.002 

N»52. Mean 

Rain at Control Station: 

0.005 inch. 

♦Inches. 


Table S-S. 

Scores for Storm September 

17-18, 1980 


Old Curve 

Q 

G 

Qa 

Qm 

Qam 

Correlation 

Coefficient 

0.133 

0.619 

0.642 

0.605 

0.637 

Mean Absolute 
Error* 

0.351 

0.227 

0,232 

0.260 

0.229 

Mean Algebraic 
Error* 

-0.307 

-0.051 

-0.060 

-0.179 

-0.055 

New Curve 






Correlation 

Coefficient 

0.276 

0.619 

0.637 

0.609 

0.631 

Mean Absolute 
Error* 

0.316 

0.227 

0.229 

0.251 

0.228 

Mean Algebraic 
Error* 

-0.246 

-0.051 

-0.057 

-0.148 

-0.054 


N«56. Mean Rain at Control Stations: 0.364 inch. 


♦Inches . 


S-S 
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Tabl* 5*6. Scorts fqr Stoni S«pttBb«r 24-2S, 1980 


Old Curve 

Q 

G 

Q« 

Qm 

Qaa 

Correlation 

Coefficient 

0.009 

0.501 

0.326 

0.298 

0.425 

Mean Absolute 
Error* 

0.397 

0.2S2 

0.280 

.. 0.287 

0.262 

Mean Algebraic 
Error* 

-0.7,93 

-0.031 

-0.027 

-0.162 

-0.029 

New Curve 






Correlation 

Coefficient 

-0.026 

0.501 

0.405 

0.384 

0.462 

Mean Absolute 
Error* 

0.323 

0.252 

0.265 

0.270 

0.258 

Mean Algebraic 
Error* 

-0.112 

-0.031 

-0.022 

-0,071 

-0.026 

N*56 . Mean 

Rain at Control 

Station : 0 

.434 inch. 

* Inches. 


Table 5-7. Scores 

for Storm 

September 25 

-26, 1980 


Old Curve 

Q 

G 

Qa 

Qm 

Qan 

Correlation 

Coefficient 

-.00002 

0.263 

0.057 

0.032 

0.152 

Mean Absolute 
Error* 

1.014 

0.071 

0.115 

0.512 

0.089 

Mean Algebraic 
Error* 

+1.013 

+0.010 

+0.037 

+0.512 

+0.024 

New Curve 






Correlation 

Coefficient 

0.015 

0.263 

0.205 

0.062 

0.252 

Mean Absolute 
Error* 

1.1S4 

0.071 

0,089 

0.532 

0.078 

Mean Albebraic 
Error* 

+1.154 

+0.010 

+0.022 

+0.582 

+0.016 


Mean Rain at Control Stations: 0.068 inch. ^Inches. 
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Tablt 

5-8. 

Scores for Storm October 24-25 

, 1980 




Q 

G 

Qa 

Qm 

Qam 

Old Curva 







CorreUtion 

Coefficient 


0.429 

0.575 

0,575 

0.556 

0.583 

Mean Absolute 


0.504 

0.268 

0,288 

0.316 

0.268 

Error* 







Mean Algebraic 
Error* 


-0.282 

+0.054 

+0,075 

-0.114 

+0.064 

New Curve 







Correlation 

Coefficient 


0.456 

0.S7S 

0.590 

0.578 

0.588 

Mean Absolute 
Error* 


0.377 

0.268 

0.275 

0.289 

0.266 

Mean Algebraic 
Error* 


-0.106 

+0.054 

+0.077 

-0.026 

+0 . 066 

N>49 . Mean 

Rain 

at Control Station: 

0.839 inch. * 

Inches . 



Table S-9. 

Scores for Storm November IS- 

16, 1980 


Old Curve 


Q 

G 

Qa 

Qm 

Qam 

Correlation 

Coefficient 


0.111 

0.530 

0.240 

0.148 

0.288 

Mean Absolute 
Error* 


0,457 

0.060 

0.304 

0.247 

0.177 

Mean Algebraic 
Error* 


+0.340 

+0.002 

0.237 

+0.171 

+0.120 

New Curve 







Correlation 

Coefficient 


0.093 

0.530 

0.225 

0.125 

0.268 

Mean Absolute 
Error* 


0.884 

0.060 

0.337 

0.446 

0.194 

Mean Algebraic 
Error* 


+0.884 

+0.002 

+0.269 

+0.443 

+0.135 

N*48. Mean 

Rain 

i at Control Stations; 

0.119 inch. 

* Inches . 
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Table 

VI 

1 

o 

e 

Scores for Stoxm November 

17-18, 1980 


Old Curve 

Q 

G 

Q» 

Qm 

Qam 

Correlation 

Coefficient 

0.363 

0.682 

0.SO6 

0.591 

0.606 

Mean Absolute 
Error* 

0.477 

0.186 

0.221 

0.280 

0.202 

Mean Algebraic 
Error* 

-0.428 

+0.018 

+0.031 

-0.205 

+0.025 

New Curve 






Correlation 

Coefficient 

0.42S 

0.682 

, 0.S51 

0.608 

0.625 

Mean Absolute 
Error* 

0.286 

0.186 

0.210 

0.221 

0.196 

Mean Algebraic 

-0.135 

+0.018 

+0.018 

-0.058 

+0.018 

Error* 






N«52. Mean Rain at 

Control 

Station: 0.764 

inch. 

* Inches . 


Table 

S-11. 

Scores for Storm 

December 

9-10, 1980 


Old Curve 

Q 

G 

Qa 

Qm 

Qam 

Correlation 

Coefficient 

-0.08i 

0.197 

0.078 

-0.056 

0.129 

Mean Absolute 
Error* 

1.663 

0.096 

0.213 

0.836 

0.136 

Mean Algebraic 
Error* 

+1.663 

+0.009 

+0.061 

+0.836 

+0.035 


New Curve 


Correlation 

Coefficient 

-0.136 

0.197 

0.006 

-0.081 

0.077 

Mean Absolu.;ie 
Error* 

1.638 

0.096 

0.170 

0.824 

0.120 

Mean Algebraic 
Error* 

+1.638 

+0.009 

+0.041 

+0.824 

+0.025 

N>53 . Mean 

Rain at Control 

Stations : 

0.212 inch. 

♦Inches. 
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Table 5-12. 

Scores for Storm March 4-5, 

, 1981 


Old Curve 

Q 

G 

Q« 

Qm 

Qam 

Correlation 

Coefficient 

-0.027 

0.726 

0.207 

0.063 

0.297 

Mean Absolute 
Error* 

1.112 

0.037 

0.494 

0.560 

0.274 

Mean Albebraic 
Error* 

♦1.112 

+0.007 

+0.315 

+0.559 

+0.161 

New Curve 






Correlation 

Coefficient 

-0.018 

0.726 

0.210 

0.066 

0.291 

Mean Absolute 
Error* 

1.205 

0.087 

0.533 

0.606 

0.289 

Mean Algebraic 

+1.205 

+0.007 

+0,362 

+0.606 

+0.184 

N"57. Mean 

Rain at Control 

Station: 

0.448 inch. 

♦Inches. 

* 


Table S-13. Scores for 

Storm March 5-6, 

1981 


Old Curve 

Q 

G 

Qa 

Qm 

Qam 

Correlation 

Coefficient 

0.257 

0.532 

0.545 

0.564 

0.539 

Mean Absolute 
Error* 

0.041 

0.046 

0.045 

0.033 

0.045 

Mean Algebraic 
Error* 

-0.038 

+0.027 

+0.025 

-0.005 

+0.026 

New Curve 






Correlation 

Coefficient 

0.030 

0.532 

0.495 

0.386 

0.524 

Mean Absolute 
Error* 

0.149 

0.046 

0.051 

0.094 

0.047 

Mean Algebraic 
Error* 

+0.134 

+0.027 

+0.026 

+0,081 

+0.026 

N«54>. Mean 

Rain at Control 

Stations ; 

0.041 inch. 

♦Inches. 



5-9 



T«bl0 S-I4 


Scores for Stom March 16-17, 1981 


Old Curve 

Q 

G 

Q» 

Qm 

Qan 

Correlation 

Coefficient 

0.0S9 

0.658 

0.482 

0.510 

0.587 

Mean Absolute 
Error* 

0.107 

0.059 

0.082 

. 0.079 

0.069 

Mean Algebraic 
Error* 

-0.004 

-0.005 

+0.005 

-0.005 

-O.QQQiS.. 

New Curve 






Correlation 

Coefficient 

0.289 

0.058 

0.604 

0,559 

0.643 

Mean Absolute 
Error* 

0.308 

0.059 

0.067 

0.164 

0.060 

Mean Algebraic 

+0,299 

-0.005 

-0.005 

+0.147 

-Q.Q05 

N*S0 . Mean 

Rain at Control 

Station : 

0.097 inch. 

* Inches. 



Table S-IS. Scores for Storm March 22-23 

, 1981 


Old Curve 

Q 

G 

Qa 

Qm 

Qam 

Correlation 

Coefficient 

-0.062 

0.805 

0.762 

0.411 

0.303 

Mean Absolute 
Error; 

1.652 

0.148 

0.182 

0.362 

0.154 

Mean Algebraic 
Error: 

+1.647 

+0.062 

+aj2SJL- 

+0.854 

+0.057 

New Curve 






Correlation 

Coefficient 

-0.022 

0.805 

0.799 

0.598 

0.309 

Mean Absolute 
Error: 

1.330 

0.148 

0.150 

0.699 

0.145 

Mean Algebraic 

+1.330 

+0.062 

+0.051 

+0.696 

+0.056 

N«53. Mean 

Rain at Control 

Stations : 

0.477 inch. 

♦Inches. 
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Table 5-16. Scores for Storm March 29-30, 1981 


Old Curve 

Q 

G 

Qa 

Qm 

Qam 

Correlation 

Coefficient 

0.432 

0.714 

0.392 

0.552 

0.616 

Mean Absolute 
Error* 

1.022 

0.113 

0.161 

0.499 

0.130 

Mean Algebraic 

+1.012 

-0.052 

+0.004 

+0.480 

-0.024 

New Curve 






Correlation 

Coefficient 

0.437 

0.714 

0.555 

0.587 

0.671 

Mean Absolute 
Error* 

1.058 

0.113 

0.137 

0.517 

0.123 

Mean Algebraic 
Error* 

N«51. Mean Rain 

+1.054 
at Control 

-0.052 

Stations: 0.196 

-0.020 

inch. 

+0.501 
* Inches . 

-0.036 


Table 5-17 

. Scores for Storm March 

30-31, 1981 


Old Curve 

Q 

G 

Qa 

Qm 

Qam 

Correlation 

Coefficient 

0.170 

0.498 

0.445 

0.456 

0.484 

Mean Absolute 
Error: 

0.475 

0.174 

0.203 

0.298 

0.184 

Mean Algebraic 
Error; 

+0.413 

-0.016 

-0.034 

+0 . 199 

-0.025 

New Curve 






Correlation 

Coefficient 

0.284 

0.498 

0.493 

0.523 

0.500 

Mean Absolute 
Error* 

0.388 

0.174 

0.182 

0.260 

0.177 

Mean Algebraic 
Error* 

+0.312 

-0.016 

-0.026 

+0.148 

-0.021 


N»S1. Mean Rain at Control Stations: 0.435 inch. * Inches. 
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Tibl* S-18. Scores for Storm April 14<1S, 1981 


1 


Old Curve 

Q 

G 

Q» 

Qm 

Qam 

Correlation 

Coefficient 

0.274 

0.S65 

0.434 

0.450 

0.543 

Mean Absolute 
Error* 

0.213 

0.112 

0.144 

0.142 

0.116 

Mean Albebraic 
Error* 

•t-0.039 

-0.012 

-0.008 

+0.014 

-0.010 

New Curve 






Correlation 

Coefficient 

0.067 

0.565 

0.448 

0.304 

0.539 

Mean Absolute 
Error* 

0.268 

0.112 

0.131 

0.176 

0.115 

Mean Algebraic 

+0.204 

-0.012 

-0.013 

+0 . 096 

-Q.Q12 

Error* 






Nb 53. Mean 

Rain at Control Station: 

0.234 inch. 

♦Inches. 



Table 5-19. 

Scores for Storm May 15-16, 

1981 


Old Curve 

Q 

G 

Qa 

Qm 

Qam 

Correlation 

Coefficient 

0.683 

0.620 

0.732 

0.741 

0.689 

Mean Absolute 
Error* 

0.067 

0.067 

0.059 

O.OSS 

0.059 

Mean Algebraic 
Error* 

-0.062 

+0.003 

-0.005 

-0.029 

-0.001 


New Curve 


Correlation 

Coefficient 

0.624 

0.620 

0.716 

0.718 

0,674 

Mean Absolute 
Error* 

0.064 

0.067 

0.057 

0.058 

0.060 

Mean Algebraic 
Error* 

-0.051 

+0.003 

-0.008 

-0.024 

-0.003 


N*53. Mean Rain at Control Stations: 0,096 inch. • Inches. 
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Table S-20. Scores for Stom May lS-19, 1981 


Old Curve 

Q 

G 

Q* 

Qm 

Qam 

Correlation 

Coefficient 

0.498 

0.544 

0.588 

0.585 

0.611 

Mean Absolute 
Error* 

0,310 

0.208 

0.213 

0.221 

0.196 

Mean Algebraic 
Error* 

4-0.049 

+0.040 

♦0.045 

+0.045 

♦0.043 

New Curve 






Correlation 

Coefficient 

0.556 

0.544 

0.625 

0.634 

0.617 

Mean Absolute 
Error* 

0.326 

0.208 

0.197 

0.237 

0.195 

Mean Algebraic 

Error* 

+0.268 

+0.040 

+0.052 

+0.154 

+0.046 

N>53. Mean 

Rain at Control Station: 

0.597 inch. 

* Inches. 



Table 5-21. 

Scores for 

Storm May 19- 

-20, 1981 



Q 

G 

Qa 

Qm 

Qam 

Old Curve 






Correlation 

Coefficient 

0.483 

0.559 

0.570 

0.628 

0.594 

Mean Absolute 
Error* 

0.143 

0.119 

0.148 

0.108 

0.119 

Mean Algebraic 

-0.063 

+0.020 

+0.067 

-0.022 

+0,043 

Error* 






New Curve 






Correlation 

Coefficient 

0.327 

0,559 

0.572 

0.585 

0.580 

Mean Absolute 
Error* 

0.182 

0.119 

0.128 

0.124 

0.118 

Mean Algebraic 
Error* 

+0.137 

+0.020 

+0.040 

+0.078 

+0.030 

N»53. Mean Rain at Control 

Stations : 

0.499 inch. 

* Inches. 
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Table S-22. 

Scores for 


Q 

G 

Old Curve 



Correlation 

Coefficient 

0.389 

0.788 

Mean Absolute 
Error* 

0.278 

0.138 

Mean Algebraic 
Error* 

+0.174 

+0.062 

New Curve 



Correlation 

Coefficient 

0.359 

0.788 

Mean Absolute 
Error* 

0.281 

0.138 

Mean Algebraic 
Error* 

+0.171 

+0.062 

N»17. Mean 

Rain at Control Station: 


Storm July 21n22, 1981 



Qm 

Qam 

0.575 

0.531 

0.695 

0.168 

. 0.183 

0.153 

+0.056 

+0.118 

+0.059 

0.486 

0.509 

0.649 

0.183 

0.135 

0.161 

+0.058 

+0.117 

+0.060 

0.158 inch 

* Inches . 
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Table 5-23. 


Best Score Fi'equencies, Storms of 1980^1981 


( 


Old Curve 

Q 

G 

Qa 

Qm 

Qam 

Highest Correlation 
Coefficient 

0 

14 

1 

3 

3 

Least Mean 
Absolute Error 

0 

16 

0 

3 

2 

Least Mean 
Algebraic Error 

New Curve 

0 

11 

6 . 

1 

3 

Highest Correlation 
Coefficient 

d 

13 

2 

4 

2 

Least Mean 
Absolute Error 

0 

15 

1 

0 

5 

Least Mean 
Algebraic 

0 

11* 

6* 

1 

3* 


★ 


Includes ties 
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Precipitation for the storm of March 4"5, 1981 was more than double that 
for April 14-15. Again the distribution was quite uniform. CSee Figure 
S-2.) The G method easily beat the satellite methods in this storm. 

The storm of October 24-25, 1980 was fairly uniform, but there were definite 
exceptions. (See Figure B>3.j Deerfield CO. 05 inch) is surrounded by stations 
with 0,54, 0.54 and 1.22 inches. Free Union (0.10 inch) is surrounded by 
stations with 0.54, 1.27, 0,91, 1.10, 1.33 and 2.08 inches. Qan verified 
best in this storm. 

The storm of August 15-16, 1980 demonstrates the irregular patterns of summer. 
(See Figure 5-4.) Thunderstorms in southeastern Virginia and the Eastern 
Shore brought more than an inch in many places, and 4.13 inches to Norfolk. 

A gage near Virginia Beach received six inches in less than two hours. 
Significant rain occurred in irregular patterns in northern and north central 
Virginia, but most of the James River Basin above Richmond got no rain. The 
highest scores in this storm were made by Qam. Thus, as the rainfall patterns 
become more irregular, the satellite's contribution to the technique increases 
significantly. 

Overall the highest correlations were obtained in the summer storm of August 
15-16, but R was also relatively high in the storms of September 9-10, 1980 
and May 15-16, 1981. The poorest correlations occurred on December 9-10, 
September 16-17, and September 25-26 (all 1980). Significantly, the mean of 
all correlation coefficients for Q, Qa, Qm and Qam for the colder period from 
September 16 to March 17 is a poor 0.300, while for the warmer period from 
March 22 to September 10 the mean is 0.560 (not that these dates are ideal 
cutoffs between warm and cold weather, but they do dramatically separate the 
lower scores from the higher.) These statements hold true regardless of which 
temperature vs precipitation curve is used. For the G method, however, 

R = 0,520 during the colder period, and ^ = 0.675 during the warmer period. 

In the colder period the G method has the highest R value in 73 percent of 
the storms when the old curve is used, and in 82 percent of the storms when 
the new curve is used. In the warmer period G has the highest R in 60 percent 
of the storms with the old curve, and only 40 percent with the new curve. 

In other words, the G method appears to be fairly reliable throughout the 
year, while the methods involving satellite data are competitive in the 
warmer months but suspect in the colder. 
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All methods overestimated rainfall in 1980 and 1981, while all methods under- 
estimated the rain in 1978. Yet a comparison of Table 2-1 and Figure 3-2 
would Suggest the opposite. That is, the cloud top temperature vs. precipi- 
tation relationship used in 1973 should produce slightly higher estimates 
than either curve used in 1980 and 1981. The answer to this apparent paradox 
is that the period March to August 1978 had above normal. precipitation while 
June 1980 to April 1981 was unusually dry. 

In regions with very dry environments (deserts and temporary drought areas) 
disappointingly scant rain is produced as a rule by the occasional cumulonimbus 
that does occur, even when the convective chimney is quite large and tall. 

In very wet environments, on the other hand, rather innocuous appearing 
build-ups may bring surprisingly heavy rain. There is usually a lag of several 
months from the onset of the below normal rainfall which initiates the drought 
condition to the inhibiting effect on such cumulonimbus as do occur. A 
similar lag occurs in the e.xcessively wet periods. 

Departures from normal in precipitation for stations representing si.x geo- 
graphical areas of Virginia are presented in Tables 5-24 and S-25. Table 5-24 
shows well above normal rainfall for the period March 1978 through August 
1978, and below normal for the following two months. Table 5-25 shows large 
deficits for the period March 1980 through April 1981. The cumulative deficit 
ranges from -2.28 inches at Wytheville to -IS. 37 inches at Norfolk. Under 
these circumstances any predictive curve ba.sed on normal precipitation will 
underestimate in the wet periods and overestimate in the dry. However, an 
examination of the scores in Tables 2-2 through 2-9 , and 5-1 through 5-22 , 
will show that Q and Qra are the only methods that are seriously affected by 
these anomalous situations. 

For the seven storms of 1978 taken together, the mean underestimate for Q is 
-0.20 inch, for Qra is -0.11 inch, and for G, Qa and Qam a negligible -0.01 
inch. The worst score for a single storm is -0.69 inch for Q, -0.38 inch for 
Qm, -0.20 inch for G, -0.19 inch for Qam, and -0.17 inch for Qa. The large 
underestimates for Q and Qm occurred on May 4-5, 1978; those for G, Qam and 
Qa on July 31-August 1, 1978. 

For the 21 storms of 1980-81 taken together, the mean overestimate for Q is 
+0.31 inch (old curve) and +0.43 inch (new curve); for Qm it is +0.16 inch 
Cold curve) and +0.22 inch Cnew curve). Qa and Qam show negligible over- 
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Tabl« 5-24. Departure from Nomal in Precipitation, 1978 (inches) 


1978 

RIC 

WY 

LYH 

ROA 

ORF 

lAD 

March 

+2.29 

+0.90 

+0.38 

+2,59 

+4.38 

+2,85 

April 

+1,54 

+1.90 

+2.60 

+4.74 

+0.19 

-1.34 

May 

+0.50 

+0.26 

+3.47 

+1.38 

+2.30 

+1.37 

June 

+1.74 

-1.45 

-0.07 

-1.46 

+4.22 

+0.75 

July 

-1,39 

-0.27 

+0.78 

+0.93 

-1.51 

+0.40 

August 

+0.87 

+4.13 . 

-0.16 

+2.18 

-4.26 

+0.04 

Cumulative 

March-August 

+5.55 

+5,47 

+7.00 

+10.36 

+5.32 

+4.07 

September 

-3,32 

-2.90 

-3.28 

-2.97 

-3.03 

-2.51 

October 

-1.73 

-1.87 

-1.76 

-2.41 

-1.56 

-1.95 

Cumulative 

September- 

October 

-5.05 

-4.77 

-5.04 

-5.38 

-4.59 

-4,46 

Legend: 

RIC * Richmond, Eastern Piedmont 

IVY » Wytheville, Southwestern Mountain 




LYH * Lynchburg, Western Piedmont 

ROA > Roanoke, Central Mountain 

ORF * Norfolk, Tidewater 

lAD ■ Dulles International Airport, Northern 
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Dtpartuns from NotimI in Pracipitatlon, 1980-1 Cinches) 


Table 5-25. 


1980 

RIC 

WY 

LYH 

ROA 

ORF 

lAD 

June 

-3.14 

-1.06 

-2.7S 

-1.70 

-2.23 

-1.72 

July 

-0.4S 

♦•2.93 

-0.44 

♦ 1.44 

-3.85 

♦0.29 

August 

-2.91 

-2.25 

-2.71 

•1.28 

'■ -1.38 

-2.58 

September 

-1.21 

♦0.27 

-1.51 

-1.76 

-2.73 

-0.59 

October 

♦■4.02 

♦0.67 

-0.25 

-0.39 

♦1.15 

♦0.03 

November 

-1.02 

-0.65 

♦0.19 

-0.70 

-0.93 

+0.36 

December 

-2.82 

-1.98 

-2.65 

-2.51 

-0.47 

-2.79 


1981 


January 

-2.22 

-1.47 

-2.28 

-2. 45 

-2.30 

-2.44 

February 

-0.25 

-0.05 

♦1.02 

-0.66 

-1.05 

+ 1.49 

March 

-1.86 

-1.14 

-1.65 

-1.03 

-1.54 

-2.49 

.i.pril 

♦0.19 

-0.36 

-0.29 

-1.05 

-0.45 

+0.09 

Cumulative 
June- April 

-11.67 

-5.09 

-13.35 

-12.59 

-15.78 

-10.35 


Legend; RIC » Richmond, Eastern Piedmont 

WY a Wyxheville, Southwestern Mountain 

LYH ■ Lynchburg, Western Piedmont 

ROA ■ Roanoke, Central Mountain 

ORF ■ Norfolk, tidewater 

lAO ■ Dulles Intemationl Airport, Northern 
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estimates of 4-0.04 and >0.02 inch for each curve i, and G has zero bias. The 
largest overestimates for Q are >1.66 inches (old curve) and >1.64 inches 
(now curve) on December 0-10, 1980, The largest for Qm are >0.85 inch (old 
curve, March 22-23, 1981) and >0,82 inch (new curve, December 9-10, 1980). 
Largest for Qa are >0,31 inch (old curve) and >0,36 inch (new curve) both on 
March 4-5, 1981. Largest for Qam are >0.16 inch (old curve), and >0.18 inch 
(new curve) both on March 4-5, 1981. G has no bias larger than >0.06 inch 
in any of the 21 storms, an excellent record. These statistics show that 
the automated method has a built-in protection against bias. 

Neither the old nor the new curve of cloud top temperature vs. precipitation 
shows decided superiority. The frequencies in Table 5-23 give a slight edge 
to the new curve. Tables S-2 through 5-22 indicate that Q was improved in 
only 38 percent of the cases when the new curve replaced the old, and Qm in 
only 33 percent of the cases. On the other hand, Qa improved in 67 percent 
of the cases, and Qam improved in 71 percent of the cases. However, the 
composite scores in Table S-lshow much larger mean errors for Q and Qm under 
the new curve than under the old, but very small iinprovements in Qa and Qam 
under the new curve. This tends to swing the verdict in favor of the old 
curve . 

Since neither curve has demonstrated marked superiority when tested in 21 
storms, it is fortunate that both curves are now integral parts of the 
computer program. The user may g.fiopt either or both curves for each rainy 
day, at his discretion. 


SECTION 6- THE FUTURE 


This report completes the reseaz'ch and development phase of the project. The 
automated technique is now ready for the operational phase for which it was 
designed. In presenting reliable precipitation animates for every cell in 
Viginia in near real time on a daily on-going basis, the technique requires 
on the order of 12S to ISO daily gage readings by dependable, highly motivated 
observers distributed as uniformly as feasible across the state. 

Appendix C lists the stations used in the various stages of the investigation. 
These include official weather stations which report daily over national 
teletype and facsimile circuits; NATS stations recruited by VPISSU, also 
reporting daily; cooperative observers whose reports are difficult to obtain 
prior to publication twomonths or more after the fact; and a handful of 
stations no longer reporting. Observation time for all these stations is, of 
necessity, between 6 a.m, and 9 a.m. local time. For best results with the 
automated technique, the cooperative observers should either be induced to 
assume NATS status or be by-passed with nearby recruits. The cooperative 
observers are in general a dedicated breed, yet their reports are of little 
value to the program if delayed more than 24 hours past observation time. 

The motivation of volunteer observers depends largely on the "human approach", 
on convincing the observer that his contribution is important, and maintaining 
that conviction by frequent feedback of positive results achieved by the 
system. The NATS observers, recruited mainly from the agricultural community, 
no doubt have a stake in the economic benefits of the program, but pride in 
having a worthwile impact on their neighbors' management problems and 
programs also looms large. 

John F. Moses Gooses has computerized the most time-consuming and laborious 
aspects of the Scofield/Oliver scheme and developed an interactive man-machine 
technique. Tests of his model are very promising; particularly good scores 
were made in a heavy thunderstorm regime over a high-density network of rain 
gages in northern Illinois. This interactive tehnique is still used almost 
exclusively for heavy rain situations, but it is hoped that eventually it can 
be applied to weaker storms, as required by agricultural users. 
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Project AGRISTARS, funded by Congress to devise a comprehensive program to 
estimate precipitation on a world-wide, year-round basis, has been examining 
all existent methods to attain this end. The Applications Laboratory of NESS, 
Washington, 0. C. is developing meteorological satellite derived products in 
support of AGRISTARS. The users of the automated technique might find it 
advantageous to maintain contact with the Applications Lab in order to benefit 
by its findings. 

In order to exploit new breakthroughs , the software of the system described 
in this report is adaptable and open-ended. The hardware, unlike many systems 
designed for short-range research projects, has been provided by NASA for a 
permanent operation, >and is generally of highest quality and durability. 
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Official Synopses by NWS Forecast Office. Washington. D. C. 

August lS-16, 1980: A cold front extension from Central New York to Indiana 

on the morning of August IS moved southeastward across Virginia in the night 
and became stationary over the Carolines on the 16th. 

September 9-10, 1980: A strong cold front located in the Ohio Valley the 

morning of September 9 moved across the Appalachians by night reaching the 
north portion of Chesapeake Bay late at night. It moved east of the Bay and 
east coast early Wednesday CSep. 10) morning, but trailed back across western 
North Carolina at 0939Z September 10. 

September 16-17, 1980: No front in the area. Southeasterly winds of 10 to 

IS knots in the afternoon becoming southerly 10 to 15 knots in the night. 

September 17-18, 1980: A cold front moved eastward through the Ohio Valley 

on September 17 and stalled over the Appalachians late in the day, spreading 
showers and thunderstorms over Virginia. A warm front formed across southern 
Virignia during the afternoon, but weakened in the night. The cold front 
extended from eastern New England through southwest Virginia by 09392 
September 18. It moved southeastward off the coast by late in the day. 

September 24-2S, 1980: A stationary front extended east-west across the. <• 

Carolinas near 35“ north. Some rain developed across Virginia due to the 
nearness of the front. 

September 25-26, 1980: A cold front extending from Ohio into Tennessee on 

the 25th moved rapidly eastward crossing Chesapeake Bay late oh the night of 
September 25, A strong secondary cold front crossed the Appalachians on 
the morning of September 26 and moved offshore early in the afternoon. 

October 1-2, 1980: A low center east oi” the Virginia capes moved slowly 

east northeastward on October 1 and early Octobea^’ 2. A cold, front in the 
midwest moved eastward across the Appalachians late on October 2. 

October 24-25, 1980: A deepening low off the Georgia coast on the morning 

of October 24 moved northward, reaching the North Carolina coast south of 
Hatteras at 3 a.m. (EDT) October 25. It continued to intensify as it moved 
northward into the southern Chesapeake Bay. 

November 14-15, 1980: A slow-moving cold front crossed the Appalachians on 

the afternoon of November 14 and into central Virginia at night. 
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November 15-16, 1980: The cold front moved across the Chesapeake Bay and 

southern Virginia the morning of November 15 and became stationary over the 
Carolines by night. Late on the night of November 15 a low moved eastward 
off the North Carolina coast. High pressure dominated Virginia on the morning 
of November 16. 

November 17-18, 1980: A frontal trough developed late the morning of November 

17 near the Virginia capes and Carolina coastline, intensifying in the after- 
noon. A low pressure center formed offshore along the front near southeastern 
Virginia during the night and moved east northeastward at 30 knots, with a 
cold front trailing down across South Carolina by late on November 18. 

December 9-10, 1980; A low over Kentucky on the morning of December 9 moved 
northeastward across West Virginia in the afternoon, across Pennsylvania in 
the night and off the New England coast on December 10. A cold front south- 
ward from the low crossed Virginia in the afternoon and night of the 9th. 

By 1200Z December 10 the front was in southeastern Virginia moving offshore. 

March 4-5, 1981; An intensifying low over the plains states was moving east- 
ward, but could not affect Virginia weather. Another intensifying low 
developed over South Carolina in the early morning of March 5 and began moving 
northeastward. 

March 5-6, 1981: The intensifying low over South Carolina at 0939Z March 5 

moved to northeastern North Carolina by 1S39Z and off the Virginia capes by 
2139Z and moved off to the east northeast thereafter, Virginia was dominated 
by strong northerly winds, with snow flurries. 

March 16-17, 1981: A strong cold front moved southeastward across Virginia 

on March 16, and out over the Atlantic late in the day. 

March 22-23, 1981: A weak high pressure ridge remained over Virginia on 

March 22, moving off the coast early in the morning of March 23. A de- 
veloping low over Arkansas early on March 22 moved rapidly eastward and off 
the South Carolina coast by late in the morning of March 23. 

March. 29-30, 1981: -A strong high off the Carolina coast moved slowly east- 
ward. A cold front over the great plains moved rapidly eastward but was 
still west of the Appalachians early on the morning of March 30. 

March 30-31, 1981; The cold front extending southward from Chicago on the 
morning of March 30 moved eastward across Virginia late in the afternoon and 
off the east coast at night. High pressure built up behind the front. 

B-3 



April 14-lS, 1981; A warm front lyings across southwest Virginia in the 
morning moved eastward across the state in the afternoon. A strong cold 
front followed in the wake of the warm front, moving east and south of 
Virginia shortly before midnight. 

May 15-16, 1981: A cold front from Ohio to Georgia in the morning moved to 

a line from eastern Ohio to central South Carolina by evening. It crossed 
Virginia during the night reaching the coast shortly after daybreak. 

May 18-19, 1981: High pressure dominated Virginia today, but a low moved into 

western Tennessee in the late afternoon with a warm front extending eastward 
from the low across the Carolinas. The front lay east-west across southern 
Virginia during the night, and became stationary across northern North Carolina 
in the early morning, 

I 

May 19-20, 1981: An east-west front across central North Carolina remained 

stationary today as a 1007-mb low in western Tennessee at 1200Z moved 
eastward along the front, reaching the North Carolina coast Wednesday morning. 
Early Wednesday morning (May 20) the front became a cold front and moved 
southward across South Carolina. 
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APPENDIX C ^ STATIONS USED IN INVESTIGATION 


Name 

Call 

ID 

Row 

Col 

Lat 

Long 

Status 

Allisonia 

AL 

1 

66 

40 

36.90 

80.75 

Control 

Altavista 

AV 

2 

63 

60 

37.10 

79.30 

Analysis 

Amelia Csee JE) 

AE 

1S3 

61 

78 

37.30 

78.03 

Analysis 

Amissville 

AM 

3 

45 

79 

38.68- 

78.02 

Control 

Appomattox 

AP 

4 

60 

66 

37.37 

78.88 

Analysis 

Bedford 

BO 

S 

60 

57 

37.35 

79.52 

Control 

Berryville 

BE 

6 

40 

79 

39.15 

77.98 

Control 

Big Meadows 

BM 

7 

47 

73 

38.52 

78.43 

Analysis 

Blacksburg 

BL 

3 

62 

44 

37.18 

80.42 

Analysis 

Bland 

BN 

9 

63 

35 

37.10 

81.10 

Control 

Bremo Bluff 

BB 

10 

56 

74 

37.70 

78.30 

Analysis 

Bridgewater 

BW 

140 

49 

65 

38.33 

73.97 

Analysis N 

Brookneal 

BR 

11 

64 

65 

37.03 

78. 9S? 

Control 

Buchanan 

BU 

12 

S3 

55 

37.53 

79.68 

Analysis 

Buckingham 

BK 

13 

58 

71 

37.55 

78.55 

Control 

Buena Vista 

BV 

14 

56 

60 

37.73 

79.35 

Control 

Byllesby 

BY 

IS 

67 

36 

36.30 

80.98 

Control 

Camp Pickett 

CP 

16 

64 

79 

37.03 

77.95 

Control 

Charlotte C. H. 

cc 

17 

64 

69 

37,07 

78.65 

Control 

Charlottesville 

CT 

18 

53 

71 

38.03 

78.52 

Control 

Chari tsvl Airpt 

CHO 

117 

52 

72 

38.13 

78.45 

Analysis 

Chase City 

CY 

141 

66 

72 

36 . 33 

78.47 

Control N 

Chatham 

CM 

19 

67 

59 

36.32 

79.40 

Control 

Clarksville 

CL 

20 

69 

71 

36.62 

78.57 

Control 

Clifton Forge 

CF 

21 

55 

S3 

37.82 

79.83 

Analysis * 

Columbia 

CU 

22 

56 

77 

37.77 

78.15 

Analysis 

Concord 

CO 

23 

61 

65 

37.28 

78.97 

Control 

Copper Hill 

CH 

24 

63 

48 

37 ,10 

30.13 

Control 

Corbin 

CN 

25 

51 

' 88 

38.20 

77.37 

Control 

Covington 

CV 

26 

55 

50 

37.80 

80.00 

Control 

Craigsville 

CR 

27 

52 

59 

38.08 

79.38 

Control 

Crozier 

C2 . 

28 

57 

82 

37.63 

77.80 

Analysis 

Cumberland 

CB 

154 

59 

75 

37.50 

78.25 

Control 

Danville Arpt. 

DAN 

36 

69 

59 

36.57 

79.33 

Analysis 
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Name 

Call 

12. 

Row 

Col 

Lat 

Long 

Status 

Danville (Bridge) 

DA 

29 

69 

59 

36.58 

79.38 

Analysis 

D.C. Nat'l Airpt* 

DCA 

95 

43 

93 

38.85 

77.03 

Analysis 

Deerfield 

DE 

30 

51 

59 

38.17 

79.37 

Control * 

Dulles Int. Arpt. 

lAD 

94 

42 

87 

38.95 

77.45 

Analysis 

Earlehurst 

EA’ 

31 

57 

47 

37.67 

80.23 

Analysis 

Edinburg 

ED 

166 

43 

70 

38.35 

78.60 

Control N 

Emporia 

EM 

32 

68 

85 

36.68 

77.55 

Control 

Farmville 

FA 

33 

61 

73 

37.33 

78.38 

Analysis 

Floyd 

FY 

34 

65 

46 

36.93 

80.30 

Control 

Fort Bel voir 

DAA 

114 

45 

90 

38.72 

77.18 

Control 

Fort Eustis 

FAF 

118 

63 

98 

37.13 

76.62 

Control 

Free Union 

FU 

35 

51 

71 

38.15 

78.57 

Analysis 

Gat bright Dam 

GT 

155 

54 

51 

37.95 

79.95 

Analysis 

Glasgow 

GG 

37 

57 

58 

37.62 

79.43 

Control 

Glen Lyn 

GL 

38 

60 

38 

37.37 

80.87 

Control # 

Gordonsville 

GD 

39 

52 

76 

38.08 

78.18 

Analysis 

Goshen 

GO 

40 . 

S3 

57 

37,38 

79.50 

Analysis 

Grundy 

GR 

41 

61 

21 

37.27 

82.08 

Control 

Hillsville 

HI 

42 

68 

40 

36.67 

80.73 

Analysis 

Holcombs Rock 

HR 

43 

59 

61 

37.50 

79.27 

Control 

Holland 

HO 

44 

68 

96 

36.68 

76.78 

Control 

Hot Springs 

HS 

45 

S3 

53 

38.00 

79.83 

Control 

Huddleston 

HU 

46 

63 

57 

37.15 

79.50 

Control 

Independence 

IN 

47 

69 

34 

36.65 

81.17 

Analysis 

Isle of IVlght 

IW 

142 

65 

97 

36.92 

76.70 

Analysis N 

Jetersville CAE) 

JE 

169 

61 

78 



N 

ilohn Flannagan Lk 

FL 

48 

62 

17 

37.23 

82.35 

Analysis 

John Kerr Dam 

KD 

49 

69 

75 

36.60 

78.28 

Analysis 

Kerrs Creek 

KC 

SO 

55 

56 

37.87 

79.57 

Control 

Keysville 

KY 

156 

63 

71 

37.17 

78.52 

Analysis 

Lafayette 

LA 

51 

62 

47 

37.23 

80.22 

Control 

Lancaster 

LN 

167 

55 

100 

37.80 

76.52 

Analysis N 

Langley AFB 

LFI 

52 

64 

102 

37.08 

76.35 

Analysis 

Leesburg 

LE 

143 

40 

85 

39,12 

77.58 

Analysis N 

Luray 

LU 

53 

45 

73 

38.67 

78.38 

Control 

Lynchburg 

LYH 

54 

61 

62 

37 . 33 

79.20 

Analysis 
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Name 

Call 

ID 

Row 

Col 

Lat 

Long 

Status 

Martinsville Pil 

MF 

1S7 

68 

52 

36,70 

79.88 

Analysis 

McDowell 

MC 

S5 

49 

57 

38.33 

79.50 

Control 

Meadows of Dan 

MD 

S6 

68 

44 

36.67 

80.45 

Control 

Millgap 

MI 

S7 

49 

54 

38.35 

79.72 

Control 

Mill Run Faim 

MR 

S3 

51 

54 

38.22 

. 79.73 

Control 

Montebello 

MB 

60 

54 

63 

37.88 

79.13 

Analysis 

Mountain Grove 

MG 

61 

52 

52 

38.10 

79.88 

Analysis 

Mt. Solon 

MS 

1S8 

49 

63 

38.35 

79.08 

Analysis 

Mount Weather 

MW 

62 

41 

80 

39.07 

77.88 

Analysis 

NAS, Norfolk 

NGU 

ISl 

65 

103 

36.93 

76.23 

Control 

New Castle 

NC 

63 

59 

49 

37.50 

80.10 

Control 

Newport 

NE 

64 

61 

43 

37.32 

80.52 

Control 

Newport 6NE 

NP 

144 

60 

44 

37.35 

80.42 

Analysis N 

Newsoms 

NS 

14S 

69 

91 

36.63 

77.12 

Analysis N 

Norfolk 

ORF 

6S 

66 

104 

36.90 

78.20 

Analysis 

North Fork Lake 

NF 

66 

63 

13 

37.13 

82.65 

Control 

North River Dam 

NR 

67 

49 

61 

38.37 

79.27 

Analysis * 

Oceania 

NTO 

116 

66 

107 

36.33 

76.02 

Control 

Orange (PE) 

PE 

71 

51 

77 

38.22 

73.12 

Control 

Palmyra 

PA 

68 

55 

75 

37.87 

78.25 

Control 

Patrick Henry 

PHF 

152 

63 

100 

37 . 13 

76.50 

Analysis 

Pedlar Dam 

PD 

69 

57 

60 

37.67 

79.28 

Analysis 

Petersburg 

PG 

159 

62 

87 

37.23 

77.40 

Control 

Philpott Dam 

PT 

70 

67 

SO 

36.78 

80.03 

Control 

Piedmont R.S.COR) 

PE 

71 

51 

77 

38.22 

78.12 

Control 

Pilot 

PI 

72 

64 

45 

37.07 

80.35 

Analysis 

Piney River 

PR 

73 

56 

65 

37.70 

79.00 

Control 

Powhatan 

PW 

160 

58 

30 

37.53 

77.92 

Control 

Quantico Marine 

.,NYG 

113 

47 

89 

38.50 

77.30 

Analysis 

Radford 

RA 

74 

63 

42 

37.13 

80.55 

Control 

Randolph 

RL 

75 

65 

69 

36.98 

78.70 

Analysis 

Rapidan 

RP 

76 

SO 

78 

38.30 

78.07 

Analysis 

Richmond 

RIC 

77 

59 

38 

37.50 

77.33 

Analysis 

Riverton 

RI 

78 

42 

76 

38.93 

78.20 

Control * 

Roanoke 

ROA 

79 

61 

51 

37.32 

79.97 

Analysis 

Rockfish 

RO 

30 

55 

68 

37.80 

78.75 

.Analysis 


Cr.4 



Name 

call 

ID 

Row 

Col 

bat 

‘ *Long 

Status 

Rocky- Mount CYa) 

RM 

81 

65 

52 

37,00 

79.90 

Analysis 

South, Boston 

SB 

161 

68 

66 

36.70 

78.88 

Analysis 

Speediirall 

SP 

82 

67 

34 

36.82 

81.17 

Control 

Staffordsville 

SF 

83 

61 

40 

37.27 

80.72 

Analysis 

Stanardsviille 

sv 

147 

so 

72 

38,27 

78.45 

Control N 

Staunton. S^P. 

ST 

84 

51 

64 

38.15 

79.03 

Control 

Sterling 

SG 

8S 

42 

86 

38.98 

77.47 

Control 

Stony Creek 

SC 

86 

66 

87 

36.90 

77.45 

Analysis 

Stuart 

SU 

162 

69 

46 

36.63 

80.27 

Control 

The Plains 

TP 

87 

43 

82 

38.90 

77.75 

Control 

Tiraberville 5E 

TI 

148 

46 

69 

38.62 

78.68 

Control N 

Trout Dale 

TD 

88 

68 

30 

36.67 

81.40 

Analysis 

Tye River 

TR 

89 

57 

65 

37.63 

78.93 

Analysis 

i^akefield 6NE 

WF 

149 

64 

94 

37.02 

76.90 

Analysis N 

Wallaceton 

WL 

90 

69 

101 

36.60 

76.43 

Analysis 

Wallops Island 

WAL 

91 

54 

115 

37.95 

75.48 

Analysis 

Warrenton 

WR 

92 

45 

82, 

38.68 

77.77 

Analysis 

Washington, Va, 

WG 

93 

45 

76 

38.72 

78.17 

Analysis 

Waverly 

WV 

IIS 

64 

92 

37.03 

77.10 

Control 

West Point 

WP 

ISO 

59 

95 

37.52 

76.83 

Control N 

Willis 

WS 

163 

66 

43 

36 , 35 

80.48 

Analysis 

Winchester 

wr 

96 

39 

76 

39.20 

78.17 

Control* 

Winterpock. 

WK 

164 

61 

84 

37.33 

77.65 

Analysis 

Woodstock 

WD 

97 

43 

71 

38.88 

78.52 

Analysis 

Wool wine 

WO 

98 

68 

46 

36.72 

80.28 

Analysis 

Wytheville 

WY 

99 

65 

35 

36.93 

81.08 

Analysis 




Out of State 




Andrews AFB, Md. 

ADW 

103 

44 

95 

38,80 

76.87 

Analysis 

Baltimore, Md. 

BAL 

100 

40 

98 

39.18 

76.67 

Analysis 

Beckley, W.Va. 

BKW 

119 

56 

34 

37.78 

81.12 

Analysis 

Bluefield, W.V. 

BLF 

106 

61 

33 

37.30 

81.22 

Analysis 

Bristol, Tenn. 

TRI 

107 

70 

16 

36.48 

82.40 

Analysis 

Charleston, W.V. 

CRW 

120 

49 

27 

38.37 

81.60 

Analysis 

Dover, Del. 

DOV 

137 

40 

115 

39.12 

75.45 

Analysis 

Elizabeth C. ,N,C, 

ECG 

108 

73 

105 

36.27 

76.18 

Analysis 

Elkins, W, Va. 

EKN 

105 

43 

52 

38.88 

79.85 

Analysis 


c-s 



Name 

Call 


Row 

Col 

Lat 

■ hsni 

Status 

Greensboro, NC 

GSO 

110 

75 

SI 

36,08 

79.95 

Analysis 

Huntington, NV 

ms 

122 

49 

14 

38,37 

32.55 

Analysis 

Matinsburg, WV 

MRB 

104 

37 

79 

39.40 

77.98 

Analysis 

Patuxent, Md. 

NHK 

101 

49 

101 

33.33 

76.42 

Analysis 

Raleigh, N.C. 

RDU 

109 

77 

68 

35.37 

78.78 

Analysis 

Rocky Mount, NC 

RNI 

112 

78 

30 

35.85 

■ 77.88 

Analysis 

Salisbury, Md. 

SBY 

102 

49 

114 

38.33 

75.52 

Analysis 

Winston-Salem NC 

INT 

111 

74 

47 

36.13 

80.23 

Analysis 


N - NATS Station 

* - Station closed during later part of this study. 

# -> Glen Lyn now takes observations at 6 p.m. 

Due to closing of stations and resulting imbalance, the status of eight stations 
was changed midway through the study. Columbia, Goshen, Independence and 
Patrick Henry were changed from Control to Analysis; Corbin, Craigsville, 

Millgap and Winchester were changed from Analysis to Control. 
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